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Abstract 

' The in-medium modifications of hadron properties are briefly discussed. We re- 

strict the discussion to the lattice QCD calculations for the hadron masses, screening 
masses, decay constants and wave functions. We review the progress made so far and 
describe how to broaden its horizon. 
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1 Introduction 

(N 

"The origin of the mass" is one of the fundamental questions in physics. The mass of one 
composite system, like nucleus, is given by the summation of the masses of its constituents 
' plus the binding energy. The proton and neutrons are composite systems too. But their 

masses can't be obtained according to this principle. We have to assume that the mass is 
mainly given by the kinetic energy of the quark and gluon fields. This energy has to be 
C<*) | smaller than the potential energy which confines the qaurks in the hadron bound state. 

It is likely that the potential energy can be modified in thermal and dense medium. This 
explanation for the origin of the mass is likely valid in the chiral limit. As soon as the 
quarks get bare masses, the chiral symmetry will be spontaneously broken. Consequently, a 
mass difference appears between hadron states with same spin but different parity quantum 
numbers. This part of the mass will disappear, when the breaking of the chiral symmetry is 
repealed at a critical temperature T c . According to the lattice QCD, the chiral symmetry 
breaking might partly be restored below T c . Therefore, the hadron masses are conjectured 
to be modified at finite temperatures and densities. 

Almost all hadron properties are extracted from the exponential decay of the two-point 
correlation functions. In the lattice QCD, the quark fields are discritized on the lattice 
sites. Any resulting operator has to conserve the lattice symmetry. Such an operator can 
be projected into the different hadron states. In present work, we restrict the discussion 
to the lattice QCD calculations. We review the progress made so far. We discuss how to 
extend this powerful tool to finite temperatures and densities. 

Understanding the modifications of the hadron properties at finite temperatures and 
densities is very essential to study the behavior of the "hadronic" matter under extreme 
conditions. It is also important in order to interpret the phenomenology of hadron pro- 
duction in heavy-ion collisions. It would not be possible to give a solid description for the 
dynamics of the phase transitions, when these modifications are not taken into account [TJ. 

The most extreme(st) modification occurs along the phase boundary. Not only a rapid 
change in the degrees of freedoms is most likely expected, but also a radical rearrangement 
of the correlation and entropy. 

2 Hadron mass spectrum on the lattice 

We aim to investigate the "origin of the mass" from first principle. As mentioned above, 
all hadron properties including the mass can be extracted from the two-point correlation 
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Fig. 1: Comparison between lattice results on light hadron spectrum at zero temperature 
and vanishing density and the experimental values 



functions. The generic meson operators are given as: 

O m (x) = r a (x)T a ^(x) (1) 

where T a p are the Dirac matrices (1, 7^, 75, 7^75, er^). The subscripts give the color 
quantum numbers. For the baryons 

O b (x) = e ijk (i>T{x)Y a piP^{x)) i> k f (2) 

Then we can write down an expression for the two-point correlation function and its 
exponential decay in the Euclidean time. For pions, the operator needed to generate the 
spectrum reads 

PHp,T) = Y, e¥£n u a @n,T) l5 d a (x n ,T) (3) 

And at time r the propagator is 

C(p,r) ee <P(p, T )pt te0)>= ^_|^ exp( _^ (p>) (4) 
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where Ah (p) represent of overlap of the operator from hadron state | h > with momentum 
p. Below, we will show that the decay constant is hidden in this parameter. Then the 
propagator can be integrated over the Euclidean time- and two spatial-directions. The 
mass will be extracted from E(p). Ef L (p) = p 2 + m| is the energy of hadron state. An 
example for the hadron mass spectrum at zero momentum is given in Fig. ^ 

To extend these calculations over the region of finite temperatures, we need temporal 
two-point correlation functions. 

< o m {x) > = -i f dii>d$du {r(x)r a ^ a (x)) (^(o)r Q ^ a (o)) t 



z . 

exp {-ipMif} - S G ) (5) 



<O m (x,r)> = , y 



\ J dU T^M-Hx^) rlpM-\0,x) exp(-5 G ) (6) 
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Fig. 2: Left panel: meson masses as function of chemical potential pa. The mass of 
pseudoscalar meson does not change with pa. The mass of vector meson drops with 
increasing pa. Right panel: the masses of baryon-type states are depicted in dependence 
on pa. No special behavior can be observed here. 



where M = 7u-D« — m q is the fermion determinant and So is the gauge action. U stands 
for the link variable on the lattice. 

The lattice calculations at finite chemical potential p q are not trivial, since the fermion 
determinant becomes complex and therefore the MC simulations are no longer applicable. 
Introducing imaginary chemical potential can solve this problem. We can also apply the 
Taylor series of the operator of interest at zero p q . 
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The first estimation for the hadron masses at finite chemical potential p is reported 
in Ref. [2]. The numerical study of two-color QCD with Wilson fermions on 4 3 x 8 lattice 
shows that the mass of the pseudoscalar meson is almost constant at finite p. The vector 
meson shows a different behavior. The mass of rho-meson decreases with increasing p. 
This behavior qualitatively agrees with the predictions from the effective models, sectional 
In the baryonic channel, there is no change observed. It is important to note that these 
results must be taken as a qualitative indication about the behavior of the masses at finite 
chemical potential. As it was not possible to define the physical units, both mass and p 
are given in units of lattice spacing a. 



3 The Debye screening mass 

In contrast to the temporal correlation functions from which we construct the hadron 
masses and decay constants, the spatial correlation functions are often analyzed on the 
lattice. The reason is the limited time direction on lattice. The screening masses can be 
constructed from the spatial correlation functions. At low temperatures, the screening 
masses are identical to the pole masses. At high temperatures, they are different. At very 
high temperatures, the screening masses are expected to approach nirT. nirT gives the 
lowest Matsubara frequency. 

The screening mass measures different quantities. It depends on the properties of the 
many-body system. Therefore, it is a useful tool to invistigate the plasma phase. The 
response of the medium to the existence of a weak perturbation, like hadron, appears 
in form of dynamical Debye screening of the long-range Coulomb potential. There are 
different definitions for the screening mass. The most common one relates it to the inverse 
Debye screening length. Also it is defined as the static limit of the longitudinal vacuum 
polarization tensor Hoo^ ~^ 0) = Mjj. 

■A 




Fig. 3: The Debye screening mass, mu(T) for 2-flavor (open circles) and quenched QCD 
(solid circles) as function of T/T c 4 j. Different mass scales are also included. 



The screening mass on lattice can be obtained from the exponential falloff of a correla- 
tor of two Polyakov loops in the spatial extension. The latter is related to the free energy 
in the limit x — > oo. 

AF(x,T) oc ^^exp[-M D (T) x] (8) 

where a s (T) = g 2 (T)/4-ir is the strong coupling. Because of the interaction with the 
medium via gluon exchange, the gluon acquires an additional mass, called chromoelectric 
mass mCE. mCE results from the finite infrared limit of the polarization tensor. It is 
proportional to g(T)T in leading order perturbation theory. 

As mentioned above, Eq. |1] can be extended to finite temperature by adding 
temperature-depending terms. The integral of the resulting correlation function over the 
Euclidean time- and two spatial-dimensions results in the screening mass. 

Efcp,T) = m 2 H +? + l[&T) (9) 

Fig. El depicts the results on the screening mass for quenched and 2-flavor QCD at zero 
chemical potential. The free energy given in Eq. |HJ can be separated into internal energy 
and entropy. One finds that the internal energy of this quark-antiquark system is related 
to the screening mass. The dependence of screening mass on the temperature Mq{T) 
measures the strength of the screening response of the medium due to inserting hadron 
in it. Different mass scales are also included. We find that the effects of the hadron mass 
on the screening response are much weaker than the effects of the energy and coupling. 



In the perturbation theory, the screening mass - next-to-leading order - reads 

m D (T) = J^> + !fM^ ln ^0L\ +CNg ^T)T + O( g HTm (10) 

where cm is a non-perturbative parameter. Therefore, it can be determined by non- 
perturbative methods, i.e. lattice. 

Including the chemical potential is a non-trivial task |2J. In leading order approxima- 
tion, the screening mass at \x q ^ takes the form 

m D {T, N ) = (Nc/3 + N f /6 + (^ q /T) 2 N f /27T 2 ) 1/2 g(T)T (11) 
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Alternatively, we can use the non-perturbative methods to calculate rri£)(T, jj, q ). We 
can use the response of the screening mass at fj, q = • As in Eq. [7J we can apply Taylor 
expansion at fixed temperature, coupling and bare quark mass to calculate mo(T, /x g ) 



m D (T, fj, q ) 



m D (T) 
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H q \ dm D (T) 
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The coefficients are determined by measuring the hadron propagators and their derivatives 
at n = by using the standard MC method. 
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A^ z is the length of the lattice in the z-direction. The results for degenerate quark masses 
arrig = 0.05 are summarized in Fig. |1J We note that the screening mass in all hadron 
channels increases with increasing the temperature. This means that the response of the 
heat path to the existence of the hadrons is enhanced as the temperature raises. Above 
T c this behavior seems to be stronger. 

The influence of chemical potential starts to set on above T c . Below T c , there is almost 
to influence. In the deconfined phase the screening mass at finite \x q is larger than its 
value at \i q = 0. The mesons are much sensitive than the baryons. It is important to note 
that we get different screening masses depending on the hadron state. This means that 
the heat bath (gluonic system) is polarized around the static charge (the hadron) and its 
parton density seems to be modified, accordingly. 




Fig. 4: The screening mass, m_o as function of T/T c . Three hadron states are shown here. 
For each hadron we depict the results at fJ, q /T = 0.0 and /J-q/T = 1.0. We set the physical 
units according to T c = 0.225 GeV. 



4 The Decay Constant 

In the continuum theory, the pion decay constant is defined as 

^/.m, = <0|u 7 475^K + (p = 0)> (13) 
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We can use the temporal correlation functions given in Eq. to calculate the decay con- 
stant for the hadron of interest. At large separations and low temperatures the pole 
contributions dominate the correlation function. 



< 0(p,r) > 



C(p) cosh(E h (p)(T - T )) 



2E h 



smh(E h (p)T ) 



(14) 



where C{p) is the residue of the current when p — > 0. Applying the GMOR relation, we 
can construct the pseudoscalar correlation function, for instance, 



< 0(t) >7r 



2 ml cosh(m 7r (r - r )) 
8m^ sinh(m 7r To) 



(15) 



There are extensive lattice calculations for the decay constants of different hadrons at 
zero temperature [T. This quantity is essential for the construction of CKM (Cabibbo- 
Kabayashi-Maskawa) unitary triangle (quark mixing matrix). Its shape is directly related 
to the properties of the weak interactions and the phenomenon of CP-violation. 

Including chemical potential requires two types of modifications. The first one deals 
with the influence of the finite density of the hadron mass. This has been discussed above. 
The second one deals with the influence on the decay constant itself. 
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The superscripts * denote the in-medium modifications. < qq >* is the quark condensates 
at finite temperature and density. In Fig. [5] we show the results from the resonance gas 
model [Sj. At least according to the quark condensates, we can easily expect that the 
decay constant has to be modified in the medium. To this end, we still need to study the 
in-medium modification of the hadron masses. 
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Fig. 5: Left panel: light and strange quark condensate as function of temperature. Right 
panel: the difference between light quark condensates normalized to the vacuum value 
calculated at a fixed baryochemical potential but different isospin chemical potentials. 
Depending on the medium, the quark condensates vanish at different critical temperatures. 



5 Wave function 

We assumed that the hadronic state given in Eqs.QEJhave to fulfill the lattice symmetry. 
But they still need a spatial dimension. The quarks have to be positioned at a certian 
distance apart from each other. Even in the confined phase, the quarks are not allowed to 



shear one location, since they are fermions. Numerically, this requirement is convenient. 
We put the quarks at different lattice sites separated by a distance r, for instance. The 
separated quarks still have to remain connected to each other. To achieve this, we have 
to add gluon links. The resulting signals for the ground state have to be dominant in the 
contrast to that for the excited states For example, the pion state can be expressed as 

\ir > = ai\qq > +a 2 \qqg > +a 3 \qqgg > H (17) 

This is equivalent to adding smeared gauge fields. The added gauge fields will build up 
a gluonic cloud around the quark-antiquark system. What would this mean, technically? 
It means that one plaquette has to be rescaled. And in the corresponding link variable 
a distance r will appear, ip(x,r). Consequently, the operators, Eqs. Q El have to be 
modified. From this smearing procedure, the correlation functions becomes oscillating in 
r. The oscillation amplitude defines the wave function of the corresponding hadron state. 

The bound states can be studied by means of Bethe-Salpeter (BS) function. The norm 

1/2 

of the BS amplitude | < O|0(O)|/i, p > | is equivalent to the parameter (p) given in 
Eq. |2 The BS amplitude gives the probability of finding a pre-defined state. This new 
(smeared) operator has to remain gauge invariant. The gauge invariant BS amplitude for 
one pion with momentum p is 



A n (p,x) = <0\d(0) l5 g(d,x)u(x)\ir(p)> (18) 

where Q (0 , p) gives the path-ordered product of the gauge links that connect the point 
with x. The normalized two-point correlators construct the amplitude 

A rf . r) _ < 0|J(0;r) 75 g(0,f;r)n(f;r) u(y; 0) l5 d(y; 0)|0 > 
<0|d(0;T) 75 «(x;r) u(y; 0) l5 d(y; 0)\0 > 

From this procedure, it is now obvious, that the wave function is also constructed 
from correlation function of two operators. The first one represents the hadron state. 
The second one gives the annihilation process. This picture is supported by the smearing 
process, as we described above. Implementing this procedure under gauge invariance 
assumption straightforwardly leads to calculating the wave function on lattice. 

The wave function is an essential quantity to investigate the bound state. The influence 
of thermal and dense medium on the wave function would reflect to what extend can the 
bound state survive. On the other hand, one needs to elaborate the bound state potential 
in such a medium. We know from perturbative and non-perturbative calculations, that 
the heavy quark potential gives a solid description for the effective potential. It has been 
found that the effective potential of two quarks with a short separation is a Coulomb-type. 
When the quarks are pulled out to a very large distance, the effective potential becomes 
linear. 



6 Heavy-ion Collisions 

In this section, we discuss mass scaling laws which explain on how the strong cou- 
pling change in dense and thermal medium. Using effective Lagrangian which assumes 
spontaneously breaking of chiral symmetry, Brown and Rho suggested a mass scaling in 
medium JIU] , in order to explain the large mass enhancement of the lepton pairs observed 
in CERES experiment. 

w I* „ ( < ^ > * 

m p U \<QQ> 
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Subsequently, the reduction of the masses of baryons and mesons are equally scaled; m* B = 
xrriB and m* M = xmy[. The validity of this scaling law is restricted within the confined 
phase, especially in the region after the freeze out. The gradual restoration of the chiral 
symmetry breaking is essential to estimate this scaling law. 

Based on QCD sum rules Hatsuda and Lee suggested another scaling law. 
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where ub is the baryon density and no is the density of normal nuclear matter, 
parameter a has a value ranging between 0.15 and 0.18. 

From an effective theory with hidden local symmetry ^2j, it is expected that the 
mass of vector meson entirely vanishes at no- At this density, the pion decay constant 
will vanish too. Furthermore, a softening of the masses of a and p mesons is expected to 
be associated with the gradual restoration of chiral symmetry breaking. 



The results given in Fig. H3 illustrate one application of these models. The invariant 
mass and width modifications in heavy-ion collisions are shown. These recent results from 
different nuclear media obviously show an access. This access can only be explained, 
when the mass is scaled according to Eq. I^Jwith a = 0.15. The first direct measurement 
of the decay channel p°(770) — > tt + it~ (branching ratio 100%) is reported in Ref. |1H| . 
Also a mass shift from p + p to peripheral A + A collisions has been observed. As a 
possible explanation for this modifications one suggests the dynamical interactions with 
the surrounding matter and the distortions of the correlations and phase space. 
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Fig. 6: The invariant mass spectrum in two different targets (media) is measured at KEK. 
The observed access can only be explained, when the modification of masses in the medium 
is taken into consideration |14j 

Another significant in-medium modification has been observed while comparing the 
production and absorption in p+p, p+A and A+A collisions. Fig.0shows the strangeness 
production and its in-medium modification. We observe an obvious modification, when 
the medium is changed. Going from p + p to A + A the medium density is significantly 
increased. It is also interesting to note that the modification depends on the hadron of 
interest. For K~, the modification is larger than for K + . The reason is the different 
influences of the medium on the particle production and absorption. The production of 
K~ mainly depends on whether there are hyperons in the medium. The absorption of 
both kaons depends on the number of nucleons. 
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Fig. 7: Particle production of K~ (top panel) and K + (bottom panel) as function of the 
excess energy [T^]. The excess energy represents the difference between y/SNN and the 
energy required for K~ {K + ) production. 



7 Summary 

We have discussed the in-medium modifications of the hadron properties. We restricted 
the discussion to the lattice QCD calculations for the hadron mass, Depye screening mass, 
decay constant and wave function. We have given a brief review of the progress made so far 
and discussed how to extend the lattice calculations to finite temperatures and densities. 

There are precise lattice calculations for the masses of low-lying hadron states at zero 
temperature. Studying the thermal effects on the hadron mass is essential to understand 
the hadron production and absorption in heavy-ion collisions. Should the hadron mass 
fulfill the description given by Hagedorn in 1960's and by the dual models in the beginning 
of 1970's, our understanding of the phase diagram and the transition at high temperatures 
and densities has to be revised (modified). The dynamics controlling such a phase transi- 
tion will be different, if the hadron resonances shall keep or lose their finite volumes (finite 
wave functions) and if their masses shall be modified or remain unchanged. We might 
need to mention here that the hadron resonance gas calculations ^2] did not assume finite 
volume and mass scaling in dense and thermal medium. 

The screening mass gives the response of the medium to the existence of a weak pertur- 
bation like the statically charged hadron. On the other hand, this quantity describes the 
properties of the medium. We get information on the parton density, entropy production 
and the long-range correlations in the medium. As we have seen, there are different lattice 
calculations using different approaches. We still need to study the medium response to 
the low-lying parity partners, like a for tt, a± for p and iV(1536) for N. Also we need to 
study the screening mass of deconfined quarks in thermal and dense plasma. 

We have also discussed the decay constant. To the auther's knowledge there is no 
lattice calculations at finite temperature. At least there is no systematic study for the 
thermal effects. Using the decay constant, we can study the restoration of the chiral sym- 
metry breaking. On the other hand, the quark condensates in vacuum |17j are obviously 
modified in the medium. Including chemical potential is not trivial, because the temporal 
correlation function is related to the decay constant and simultaneously to the hadron 
mass. Both quantities have to be expressed as function of chemical potential. 

The wave function provides essential information on the bound state. There is a 
framework to proceed calculations for the wave function on lattice. We have mentioned 



Q 



that the wave function can be expressed by the Bethe-Salpeter amplitude. Studying 
the hadronic bound states at finite temperatures and densities elaborates the possible 
interactions with the medium and to what extend can the bound state survive. 

Finally we have discussed some phenomenological examples on the modification of 
some hadron properties. 
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